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Highly ionised plasma in the Large Magellanic Cloud: Evidence for 
outflows and a possible galactic wind* 
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ABSTRACT 

Based on an analysis of the interstellar highly ionised species C IV, Si IV, N V, and O VI 
observed in the FUSE and HST/STIS E140M spectra of four hot stars in the Large Magellanic 
Cloud (LMC), we find evidence for a hot LMC halo fed by energetic outflows from the LMC 
disk and even possibly an LMC galactic wind. Signatures for such outflows are the interme- 
diate and high-velocity components (i?lsr ^ 100 kms" 1 ) relative to the LMC disk observed 
in the high- and low-ion absorption profiles. The stellar environments produce strong, narrow 
(T < 2 x 10 4 K) components of C IV and Si IV associated with the LMC disk; in partic- 
ular they are likely signatures of H II regions and expanding shells. Broad components are 
observed in the profiles of C IV, Si IV, and O VI with their widths implying hot, collision- 
ally ionised gas at temperatures of a few times 10 5 K. There is a striking similarity in the 
O Vl/C IV ratios for the broad LMC and high-velocity components, suggesting much of the 
material at «lsr ^ 100 km s _1 is associated with the LMC. The velocity of the high-velocity 
component is large enough to escape altogether the LMC, polluting the intergalactic space 
between the LMC and the Milky Way. The observed high-ion ratios of the broad LMC and 
high-velocity components are consistent with those produced in conductive interfaces; such 
models are also favored by the apparent kinematically coupling between the high and the 
weakly ionised species. 
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1 INTRODUCTION 

The structure and energetics of the interstellar medium (ISM) in 
galaxies are strongly influenced by stellar winds and supernovae 
(e.g., Abbott 1982). This is particularly true for spiral and irreg- 
ular galaxies with on-going massive star formation, such as the 
Milky Way, and the Large and Small Magellanic Clouds (LMC, 
SMC). The input of energy from winds and supernovae is ex- 
pected to affect the ISM on kiloparsec scales, creating regions of 
super-heated gas (e.g., Norman & Ikeuchi 1989). The 10 6 K gas 
that expands into the halo of a disk galaxy may eventually cool 
and fall back onto the disk in a massive system, participating in a 
galactic fountain (Shapiro & Field 1976); in a low mass system, 
this material may escape the galaxy altog ether in a galactic wind 
jBreitschwerdt. Voelk. & McKenzielll99ll) . Galactic fountains and 
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winds affect the scale over which newly produced metals are dis- 
tributed in the ISM of a galaxy and/or in the intergalactic medium. 
For example, outflo ws can affect the [a/Fe] ratios inside and out- 
side dwarf galaxies dRecchi. Matteucci, & D'Ercolel200ll) . Gas dy- 
namics a re therefore criti cal in determining the evolution of galax- 
ies (e.g.. lMatteuccill2003r) . Characterising the phenomena of infall 
and outflow in nearby galaxies is crucial for understanding the role 
they may play in galaxy evolution and in the enrichment and evo- 
lution of the intergalactic medium. 

The LMC is the nearest gas-rich disk galaxy to our own Milky 
Way. As such it has been the target of intense studies of its inter- 
stellar medium (ISM), including its dust (e.g., Meixner et al. 2006), 
its molecular gas (Tumlinson et al. 2002), its warm H I (Kim et 
al. 2003), its warm ionised gas traced by Ha emission (e.g., Dan- 
forth et al. 2002), its warm-hot gas traced by O VI (Howk et al. 
2002; Danforth & Blair 2006), and its hot gas traced by X-rays 
(e.g., Snowden & Petre 1994). In this work, we report observa- 
tions of interstellar absorption from the Li-like ions N V, C IV, 
Si IV, and O VI towards stars within the LMC. These ions give us 
the means to study the highly-ionised material produced by feed- 
back processes in the Milky Way and in the Magellanic Clouds. 
Si IV, C IV, N V, and O VI peak in abundance at (0.6,1.0,1.8, 
and 2.8) x 10° K, respectively, in collisional ionisation equilibrium 
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(CIE; Sutheiiand& Dopita 1993). They are, however, unlikely to 
be in equilibrium since these temperatures correspond to the peak 
of the interstellar cooling curve, and gas at those temperatures is 
likely to cool faster than it recombines. Instead, these ions may 
trace interfaces between very hot gas (T > 10 6 K) and warm or 
cool (T < 10 4 K) material and serve as a probe of the interaction 
between the phases of the ISM having the most mass (warm/cool 
gas) and the most direct connection to feedback processes (hot gas). 

Observational searches for a hot corona about the LMC us- 
ing the highly ionised species started with WE spectra of C IV, 
Si IV , and N V towards the hot star Sk-67° 104 Jde Boer & Savage! 
1980). However, the crude spectral resolution of WE and the fact 
that that the hot st ar could be part ly responsible for producing C IV 
and Si IV brought lchu et al.l i 1994b to conclude these observations 
were not sufficient to prove the existence of a hot LMC halo (N V 
was only marginally detected at best). These highly ionised species 
could be either produced by p hotoionisation from the star itself or 
in the superbubble via shocks. Wakk er et al. | |l99l) found C IV ab- 
sorption in the spectra of cooler stars (09.7-B1) using modest res- 
olution HST spectroscopy. This was the first evidence of C IV not 
related directly to local effects, as these stars were deemed unable 
to produce significant high-ion absorption via photoionisation. 

More recently, FUSE has contributed to the effort of examin- 
ing the hot g as in the LMC by probing diffuse O VI in both M agel- 
lanic Clouds dHowk et alj20 02; Hoopes et al. 20021; iDanforth et al.l 
l2002l ; lDanforfh & Blaij|200^) . lHowk et alj J2002h used FUSE ob- 
servations of sight lines to the LMC to study the O VI con- 
tent and kinematics of the LMC. Their survey of 12 sight lines 
through the near side of the LMC probed collisionally-ionised gas 
at T ~ 3 x 10 5 K in the LMC and shows in particular that the LMC 
contains a significant amount of highly-ionised gas, with O VI col- 
umn densities equal or greater than those found along extended 
paths through the Milky Way; and much of the O VI in the LMC 
must reside in an extended "halo" or thick disk distrib ution, likely 
similar to that found in our Galaxy dSavage et alj2003h . 

These studies of the individual highly ionised species have 
broadly sketched a picture of the hot gas in the LMC and the pres- 
ence of coronal hot gas, but a crucial step is to combine all the 
available highly ionised species observed with HST/STIS E140M 
and FUSE for a better understanding of the highly ionised gas in 
the LMC. In this paper, we take such a step by presenting the de- 
tailed kinematics and high-ion absorption (O VI, C IV, Si IV, and 
N V) towards four stars in the LMC using FUSE and STIS. Two 
of these sight lines probe prominent superbubbles, while the other 
two probe H II regions, allowing us to compare the highly ionised 
species and their connection with weaker ions towards different 
types of interstellar regions, which is useful for differentiating local 
from global effects. The organisation of this paper is as follows. Af- 
ter describing the observations, data reduction, and analysis of the 
data (measurements of the kinematics, column densities, and col- 
umn density ratios) in Sj2] we present an overview of the kinematics 
in Sj3] In Sj4]we discuss the origin of the highly ionised plasma in the 
LMC and the implications of our results. A summary of the main 
results is presented in Sj5] 



2 OBSERVATIONS AND ANALYSIS 
2.1 Database 

From the data archives of FUSE and HST/STIS E140M, we se- 
lected four LMC targets. The selection criteria were to have the 



four critical highly ionised species (Si IV, C IV, N V, and O VI) 
observable in the FUV and UV and to be able to model these ions 
without too much uncertainties (i.e. with a continuum near the rel- 
evant ions that can be modeled without too much uncertainty). This 
resulted in the sample of the four stars presented in TableQ] Fortu- 
itously, these stars are also located in different interstellar environ- 
ments (from Chu et al. 1994 and references therein): 

(1) Sk— 65°22 is projected within the su pergiant shell LMC 1, 
which is about 1 kpc across dMeaburrj|"l 980). It is located in a faint 
H II region DEM 48 at the southeast rim. This supergiant shell con- 
tains the OB association LH 15. 

(2) Sk-67°21 1 is in the H II region DEM 241 that contains the OB 
association LH 82. The bright UV continuum and the early spec- 
tral type (02III(f*), Walborn et al. 2002) imply the star can likely 
produce strong C IV and Si IV absorption associated with its stellar 
environment. 

(3) Sk— 69°246 is a WR star situated in the X-ray bright superbub- 
ble 30 Dor, and more precisely is in the south end of Shell 5, which 
is the X-ray-bright "chimney" to the north of the LH 100 OB asso- 
ciation (Y.-H. Chu 2006, priv. comm.). 

(4) Sk— 71°45 is located in the superbubble DEM 221 associated 
with the OB association LH 69. 

2.2 The FUSE observations 

The FUSE observations are from the PI team programs P103 (PI: 
Sembach) and PI 17 (PI: Hutchings). Descriptions of the FUSE in- 
strument design and inflight performance are found in Moos et al. 
(2000) and Sahnow et al. (2000). The LiF 1A, LiF 2B, SiC 1A, and 
SiC 2B segment spectra cover the wavelength region of the O VI 
AA1031.926, 1037.617 doublet. However, we have chosen not to 
use the SiC 2B data because of the strong fixed pattern noise and 
relatively low resolution of observations near O VI when using this 
segment. Observations in the SiC 1A and LiF2B segments have 
~30% and ~60% of the sensitivity of the LiF 1A segment, respec- 
tively. Therefore, we principally use the LiF 1A segment for our 
measurements. We did not find any difference between the O VI 
profiles observed in the different segments (within the signal-to- 
noise). We did not combine the measurements from the three dif- 
ferent segments to maintain the optimal spectral resolution. The 
observations were obtained in the LWRS (30" x 30") apertures. We 
use the current calibration pipeline software (CALFUSE version 
3.0.8 and higher) to extract and calibrate the observations. 

The zero point in the final wavelength scale was established 
by first estimating an average FUSE velocity using the ISM lines of 
Si II A1020.699, Ar I A1066.660, Fe II A1055.262. We used these 
lines because (i) they are not too saturated to be able to measure 
the centroids accurately; (ii) they are present in the LiF 1 A segment 
and encompass the O VI A 103 1.926 interstellar absorption line, 
minimising spurious wavelength shifts; (iii) they trace the same gas 
that the ions used in the STIS wavelength range (S II, Fe II, and 
Ni II). We then shift the FUSE wavelength to the absolute LSR 
frame determined with the accurate STIS spectra. We estimate that 
the relative FUSE velocity calibration of the O VI is accurate to ~5 
km s _1 after comparing the shifts needed for the velocity centroids 
C II, Fe II, Si II in FUSE to match the LSR velocity of similar ions 
in STIS E140M. 



2.3 The STIS E140M observations 

The STIS E140M Observations were obtained through the GO pro- 
gram 8662 in Cycle 9 (PI: J. Lauroesch). This program was built 
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Table l. Summary of stellar and sightline properties 



Name 


Alias 


Type 


I 


b 


V 




u(Ha) 


Environments 








(°) 


(°) 


(mag) 


(kms -1 ) 


(kms -1 ) 




Sk-65°22 


HD 270952 


06Iaf+ 


276.40 


-35.75 


12.07 


1337 


270 ± 15 


Supergiant shell DEM 48 (N 13) 


Sk-67°211 


HD 269810 


02III(f*) 


277.70 


-32.16 


12.28 


3593 


285 ± 17 


H II region DEM 241 (N 59A) 


Sk-69°246 


HD 38282 


WN6h 


279.38 


-31.66 


11.15 


1835 


252 


X-ray bright superbubble 30 Dor 


Sk-71°45 


HD 269676 


04-5 m(f) 


281.86 


-32.02 


11.51 


2488 


227 ± 19 


Superbubble DEM 221 (N206) 



Note: Spectral types and terminal velo cities are from IWalborn et all 120021) and iMassa et all 120031) for Sk-65°22, Sk-67°211, and 
Sk- 71°45 and from [Willis et al J fe004l) for Sk-69°246. ti(Ha) are from lDanforml 120031). excep t towar ds Sk-69°246 where it comes 
from lChu et alj Jl994l) and re ferences therein. Velocities are in the LSR frame, a : from lChu et al ] 1 19941) . The designations DEM and N 
refer to entry in the catalog o f lDavies, Elliott. & Meaburi] ll976l) and lHenizd Jl956t) . respectively. 



as a snapshot survey to complement the FUSE O VI observations 
by providing coverage of the other highly ionised species, C IV, 
N V, and Si IV. The STIS data were obtained with the E140M 
intermediate-resolution echelle grating. The entrance slit was set to 
0'.'2 x 0'.'2. The spectral resolution is 7 km s -1 with a detector pixel 
size of 3.5 km s -1 . Information about STIS and its in-flight perfor- 
mance is given by Woodgate et al. (1998), Kimble et al. (1998), and 
in the STIS Instrument Handbook (Proffitt et al. 2002). 

Standard calibration and extraction procedures were employed 
using the CALSTIS routine (version 2.2). The STIS data reductions 
provide an excellent wavelength calibration with a velocity uncer- 
tainty of ~1 kms -1 . The heliocentric velocity was corrected to the 
dynamical local standard of rest-frame (LSR) (Mihalas & Binney 
1981). 

2.4 Selection of species and velocity structures 

The main goal of this work is to study the properties and ori- 
gin^) of the highly ionised species in the LMC. O VI A1031.926, 
N V AA1238.821, 1242.804, C IV AA1548. 195, 1550.770, and Si IV 
AA1393. 755, 1402.770 are the highly ionised species at hand. The 
excitation potentials are 113.9 eV, 77.5 eV, 47.9 eV, and 33.5 eV 
for O VI, N V, C IV, and Si IV, respectively. Therefore, O VI and 
N V are better diagnostics of shock-heated gas than C IV and Si IV, 
the latter being susceptible to photoionisation. 

To understand the connection, if any, between the highly 
ionised species and the weakly ionised species, we have to investi- 
gate the kinematics of the ions tracing cool/warm and hot gas. For 
tracing the cool/warm gas we use the low-ionisation species O I 
A1302.169 Si II A1526.707, S II A1250.578, and Fe III Al 122.524. 
O I A1302.169 and Si II A1526.707 are strong transitions revealing 
the weaker absorption components along each sight line, while S II 
A 1250.578 is a weaker transition that reveals the velocity structure 
of the stronger absorption components. Si II and S II trace mostly 
the warm neutral gas (WNM), although a fraction may arise in the 
warm ionised medium (WIM). O I is one of the best tracers of H I 
because the ionisation fraction of O is coupled to that of H via res- 
onant charge-exchange reactions; hence, O I should trace only neu- 
tral gas along the line of sight. Fe III Al 122.524 is found mostly in 
the WIM, and therefore is an excellent substitute for Ha emission 
measurements. 

In Figs. [T] to [4] we show the observed flux distribution in the 
vicinity of the highly ionised species with our adopted continuum 
fits (see the next sections for a full description of the continuum 
placement for those species), while in Figs. [5] and [6] we show the 
normalised profiles of the weakly and highly ionised species. The 
continuum near Si II is modeled with a straight line. The continuum 
near O I and S II are fitted with low-order Legendre polynomials. 



For the weakly ionised species, the most complicated continuum 
placement is near Fe III because of very strong stellar wind features 
from Si IV and P V. Near this transition, very high-order Legendre 
polynomials (d > 5) must be used, and this results in a large un- 
certainty in the continuum placement (see Figs.[5]and|6l. However, 
the velocity centroids of the LMC Fe III components should not be 
much affected by the continuum placement, at least near the core 
of the absorption line. 

This set of lines of sight provides an opportunity to com- 
pare the highly ionised species in various environments from the 
Galactic halo to the LMC. Insight on the highly ionised species at 
lower velocity is important in its own right, but similarities (or lack 
thereof) in the high-ion kinematics and ratios in the MW and LMC 
can provide further insight on the origin(s) of the highly ionised 
species in the LMC. In Figs.[5]and[6] absorption for the various ions 
is observed from the MW to the LMC velocities. These profiles can 
be separated naturally in four main components: (1) the MW com- 
ponent (roughly «lsr i$ 30 kms -1 ), (2) the intermediate-velocity 
component (IVC with roughly 30 < ulsr i$ 100 kms -1 ), (3) the 
high-velocity component (HVC with roughly 100 < «lsr ^ 175 
kms ), and (4) the LMC components (roughly ulsr > 175 
kms - ). Here, we follow the definition of IVC and HVC from 
Wakker (1991). We will argue that the IVC is closely connected 
to the MW, but the high- velocity component may be an HVC with 
respect to the MW or may have an LMC origin. 

2.5 Measurements of column densities, velocities, and line 
widths of the highly ionised species 

To estimate the column densities in this work, we use the atomic 
parameters compiled by Morton (2003). 

2.5.1 O VI measurements 

Our O VI measurements are based on an analysis of the 1031.926 
A transition. The O VI A1037.617 line is always contaminated by 
C II, C II*, and H2 lines, and therefore cannot be used. The O VI 
A 103 1.926 absorption in the Galactic component can be contam- 
inated by molecular lines from the Milky Way and LMC and in 
principle by CI I A1031.507. However, CI I is at -122 kms -1 and 
is only strong if the molecular fraction is high. We searched for 
the presence of CI I A 1004.668 in our four sight lines, which has a 
strength similar to CI I A 103 1.507, and found none at Galactic or 
LMC velocities. 

Three molecular lines can contaminate O VI A 103 1.926. The 
HD 6-0 R(0) A 103 1.9 12 is at -4 kms" 1 relative to O VI in the 
restframe. No absorption from HD 3-0 R(0) through 8-0 R(0) 
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Table 2. Summary of AOD measurements 



Table 3. C IV and Si IV fit results 



Ion 



(km s 1 ) 



ba 

(kms -1 ) 



log N a 
(cm" 2 ) 



Av 

(kms" 1 ) 



O VI 



N V 



258.9 ± 1.7 
128.1 ± 1.0 
53.0 ±0.5 
-3.0 ±0.6 



Sk-65' 
49.6 ± 1.6 

37.8 ± 1.2 

22.9 ± 0.6 
25.9 ± 0.4 



22 

14.14 ±0.02 
13.98 ±0.02 
13.95 ±0.02 
13.85 ±0.03 

< 13.26 

< 13.26 

< 13.20 

< 13.27 



LMC [187,350] 
HVC [88,187] 
IVC [23,88] 
MW [-45,23] 
LMC [187,350] 
HVC [88,187] 
IVC [23,88] 
MW M5.23] 



O VI 



N V 



251.2 ± 2.3 
129.1 ± 1.6 
58.2 ± 1.1 
3.8 ± 1.5 



Sk-67 c 
53.9 ± 2.6 
32.5 ± 1.3 
26.9 ±0.8 
22.2 ± 2.0 



211 

14.20 ±0.02 
13.87 ±0.02 
14.13 ±0.02 
13.97 ±0.03 

< 13.26 

< 13.12 

< 13.11 

< 13.17 



LMC [180,355] 
HVC [100,182] 
IVC [30,100] 
MW M3.30] 
LMC [180,355] 
HVC [100,182] 
IVC [30,100] 
MW M3,30] 



Sk-69°246 



O VI 



N V 



215.7 ±2.1 
103.1 ± 2.7 
62.7 ±2.2 
5.5 ± 1.4 



46.3 ± 3.5 

15.1 ±4.9 

17.4 ±2.4 

32.2 ± 2.1 



14.36 ±0.09 
13.55 ±0.20 
13.59 ±0.10 
13.96 ±0.10 

< 13.40 

< 13.20 

< 13.26 

< 13.61 



LMC [139,302] 
HVC [85,139] 
IVC [41,85] 
MW [-61,41] 
LMC [139,302] 
HVC [85,139] 
IVC [41,85] 
MW [-61,41] 



O VI 



N V 



227.3 ± 1.2 
143.0 ±0.9 
67.7 ±0.8 
-0.2 ±1.4 



Sk-71°45 
53.4 ±1.5 14.37 ±0.02 



32.5 ± 0.7 
31.2 ± 0.6 
26.0 ± 1.5 



14.12 ±0.02 
14.20 ±0.02 
13.77 ±0.02 

< 13.51 

< 13.16 

< 13.21 

< 13.22 



LMC [162,320] 
HVC [104,183] 
IVC [26,104] 
MW [-51,26] 
LMC [162,320] 
HVC [104,183] 
IVC [26,104] 
MW [-51,26] 



Note: Refer to 5 32.5. H2.5.2l for more details on these measurements. 

at 1066.271, 1054.433, 1042.847, 1021.546, and 1011.457 A was 
found in the spectrum of our target stars. Therefore, HD 6-0 R(0) 
A1031.912 does not contaminate the O VI transition. The fh lines 
(6-0) P(3) and R(4) lines at 1031.191 and 1032.356 A can also con- 
taminate the O VI absorption. These are found at —214 and +125 
km s _1 relative to the O VI restframe, respectively, and can arise in 
both the Milky Way and LMC (see Howk et al. 2002). The imprints 
of these lines can be clearly seen in the spectrum of Sk— 65°22, 
Sk-69°246, Sk-71°45, but not in the spectrum of Sk-67°211. 
To estimate the contribution to the O VI profiles from the 6-0 P(3) 
and R(4) lines of molecular hydrogen, we fitted Gaussian profiles 
to a number of J = 3 and 4 transitions of H2 with line strengths 
A/ similar to the contaminating 6-0 line transitions. The available 
H2 transitions can be found in Table 3 of Howk et al. (2002). We 
have used only H2 lines present in LiF 1A to minimise any effects 
of changes in the shape of the instrumental line spread function. 
While saturation effects may be significant, the choice of lines with 
similar A/ circumvent in great part such problems (see Howk et al. 
2002, Wakker et al. 2003). 

A difficulty in assessing the interstellar O VI absorption is that 



Ion 


V 


b 


log TV 


Cloud 




(kms -1 ) 


(km s — 1 ) 


(cm~ 2 ) 








Sk-65°22 




C IV 


271.7 ± 1.7 


23.9 ± 2.7 


13.59 ± 0.05 


LMC 




265.4 ±0.5 


6.1 ± 1.2 


13.45 ± 0.07 


LMC (H II) 




105.1 ± 13.9 


24.4 ± 18.2 


13.11 ± 0.09 


HVC 




43.1 ±0.6 


17.3 ± 0.9 


13.84 ±0.02 


IVC 




-11.0 ± 1.1 


13.2 ± 1.8 


13.35 ± 0.04 


MW 


Si IV 


274.1 ± 1.0 


24.0 ± 1.5 


13.20 ±0.03 


LMC 




269.0 ±0.3 


6.7 ± 0.5 


13.25 ±0.03 


LMC (H II) 




135.6 ±5.4 


41.4 ± 10.2 


12.59 ±0.07 


HVC 




47.3 ± 0.8 


17.9 ± 1.2 


12.96 ±0.03 


IVC 




-10.3 ±0.4 


10.7 ±0.8 


12.96 ±0.02 


MW 






Sk-67 c 


'211 




C IV 


272.2 ±0.3 


11. 4± 1.0 


(>)15.94 ± 0.32 


LMC (H II) 




251.2 ± 3.4 


43.1 ± 2.9 : 


13.80 ± 0.05 : 


LMC 




113.2 ±4.9 


41.0 ±8.3 


13.48 ±0.07 


HVC 




47.9 ±6.7 


22.9 ±6.8 


13.78 ±0.22 


IVC 




14.0 ±5.9 


23.7 ± 4.3 


13.89 ±0.16 


MW 


Si IV 


276.2 ±0.3 


10.6 ± 0.9 


(>)15.09 ± 0.25 


LMC (H II) 




262.0 ± 1.5 


37.4 ± 1.5 


13.53 ±0.03 


LMC 




38.7 ± 1.6 


31.7 ±2.3 


13.29 ±0.03 


IVC 




0.1 ± 0.7 


10.5 ± 1.2 


12.90 ±0.05 


MW 






Sk-69°246 




C IV 


214.1 ± 0.5 


16.7 ± 1.4 


(>)14.28±0.02 


LMC (H II) 




213.5 ± 1.8 


37.0 ±7.1 : 


13.90 ±0.14 : 


LMC 




98.6 ±3.5 


20.7 ± 7.9 


13.09 ±0.09 


HVC 




50.2 ± 1.3 


14.5 ± 2.5 


13.44 ±0.05 


IVC 




6.7 ± 1.0 


19.6 ± 1.7 


13.75 ±0.03 


MW 


Si IV 


218.0 ±0.3 


12.3 ± 1.2 


(>)14.23 ± 0.08 


LMC (H II) 




214.5 ±0.6 


37.0 ±7.1 


13.72 ±0.04 


LMC 




102.4 ± 2.0 


17.6 ±3.9 


12.52 ±0.06 


HVC 




51.3 ± 1.0 


16.8 ± 1.9 


12.85 ±0.03 


IVC 




3.6 ±0.5 


16.5 ± 0.8 


13.21 ±0.02 


MW 






Sk-71 


°45 




C IV 


221.0 ±0.3 


15.4 ±0.5 


(>)14.56 ± 0.03 


LMC (H II) 




128.6 : 


25.0 : 


13.35 : 


HVC 




61.4 ±0.9 


27.6 ± 1.5 


13.98 ±0.02 


IVC 




1.0 ± 1.0 


13.8 ± 1.7 


13.47 ±0.04 


MW 


Si IV 


221.2 ±0.4 


14.8 ± 0.7 


(>)14.18 ± 0.05 


LMC (H II) 




128.6 : 


25.0 : 


12.85 ±0.08 


HVC 




66.3 ± 1.2 


19.8 ±2.0 


13.34 ±0.03 


IVC 




-2.6 ± 1.2 


12.7 ±2.1 


13.05 ±0.05 


MW 



Note: The ":" symbol highlights the uncertainty in the profile fitting to the 
component where colon is marked. The "(>)" symbol indicates that the 
line is possibly uncertain due to the effects of strong saturation. See i]2.5.3| 
for more details on these measurements. 



the O VI stellar component predominantly associated with the wind 
of the star can produce a complicated continuum. Furthermore stel- 
lar wind variability can change the shape of the continuum and pro- 
duce discrete absorption components (DACs) that can mimic O VI 
interstellar absorption (Lehner et al. 2001). The matter is compli- 
cated by the fact that we cannot use interstellar O VI A1037.617 
(see above). To evaluate the effects of the stellar wind, we con- 
sidered the spectra of O VI AA1031.926, 1037.617 over a large 
velocity range in the O VI vicinity following Lehner et al. (2001). 
We note that there is no observed variation between the different 
exposures that were used to produce the final spectra. We review 
the continuum placement and possible contamination from H2 and 
stellar wind features for each star: 
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Sk~65°22: Its stellar wind has a terminal velocity Voo = 
— 1350 kms -1 , which is less in absolute value than the 1650 
kms" 1 velocity separation between the lines of the O VI doublet. 
The interstellar measurements are therefore not strongly affected 
by DACs (Lehner et al. 2001). We note that the stellar O VI profiles 
can be observed predominantly in both lines near —1000 kms" 1 . 
We are therefore confident that the stellar continuum near O VI 
A1031.926 can be simply modeled as shown in the top panel of 
Fig.Q] We also show in this figure the H2 model, where the Galac- 
tic component is weaker than the LMC component. 

Sk—67°21I: This 02 supergiant star has the simplest stellar 
continuum near O VI in our sample displaying a very prominent 
P-Cygni profile with v x = —3600 kms" 1 (see Fig.|2j. The effect 
of DACs on the local continuum is negligible since DACs are only 
strong and narrow near Voo (Lehner et al. 2001). The H2 contami- 
nation is also insignificant towards this star. 

Sk—69° 246: This is certainly the most complicated star in our 
sample for estimating the interstellar O VI. At a terminal velocity 
of Woo = -1847 kms" 1 , DACs from the O VI A1037.617 may be 
present in the galactic and LMC components of O VI A1031.926. 
Moreover, at such terminal velocity, the wind features from O VI 
A 103 1 .926 are lost in H I Lyman/9, making it impossible to identify 
features that may be present in the weaker member of the doublet 
using the stronger member. Therefore, it may be possible that the 
O VI A1031.926 interstellar profile is contaminated. We make the 
assumption it is not. The continuum placement is also uncertain. 
In Fig. [3] (top-panel) we show two likely continua. Using the H2 
modeling, the continuum cannot be lower than the solid contin- 
uum, otherwise we would underestimate the amount of the H2 in 
the LMC component. 

Sk—71°45: This star has a large terminal velocity (voo = 
—2500 kms -1 ), and therefore the interstellar O VI profile is un- 
likely to be contaminated by a DAC. The continuum is well mod- 
eled with low-order polynomial function (see Fig.Q. The Galactic 
H2 is stronger than the LMC component towards this star. 

The O VI absorption does not present distinct components, 
as can be observed in the profiles of the other ions. We adopted 
the apparent optical depth (AOD) method of Savage & Sem- 
bach (1991) to estimate the column density (N a ), velocities (v a ), 
and line widths (b a ) of O VI. The absorption profiles were con- 
verted into apparent column densities per unit velocity N a (v) = 
3.768 x 10 14 ln[F c /F obs (u)]/(/A) cm" 2 (kms" 1 )" 1 , where F c 
is the continuum flux, F ^ s (v) is the observed flux as a func- 
tion of velocity, / is the oscillator strength of the absorption 
and A is in A. The values of v a , b a , \ogN a are obtained from 
v a = JvN a (v)dv/N a ,b a = [2j(v-v a ) 2 N a (v)dv/N a ] 1 / 2 ,md 
N a = J N a (v)dv (see Sembach & Savage 1992). 

The O VI absorption profiles do not contain easily distin- 
guished absorbing components, gas associated with the LMC and 
gas at lower velocities is usually easily separated. An important 
aspect of our investigation involves studying the kinematical rela- 
tionships between the high and low ions and the column density 
ratios of the highly ionised species. We used the C IV and Si IV 
profiles to determine velocity integration limits for O VI that would 
allow us to crudely separate the different components (MW, IVC, 
HVC, and LMC). 

The measurements for O VI are presented in Table [2] For 
Sk— 69°246, we adopt the mean results of the high and low con- 
tinua. The listed ±lcr errors are from a quadrature addition of the 
statistical error and the continuum placement error. The la errors 
may be slightly underestimated (by no more than 0.05 dex in the 
log of the column density although usually lower) because they do 
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Figure 1. Plots of several highly ionised species in flux unit against the 
LSR velocity towards Sk— 65°22. The solid line is the continuum to the 
interstellar line. In the O VI panel, we also the model to the H2 line (6-0) 
P(3) and R(4) lines at 1031.191 and 1032.356 A. 

not include systematic from the removal of the H2 contamination. 
Finally, note that saturation is unlikely to affect our measurements 
since the broadening of O VI is greater than the instrumental broad- 
ening of FUSE. 

2.5.2 N V measurements 

The interstellar N V measurements are the most straightforward be- 
cause the strong stellar wind makes the continuum placement sim- 
ple. Mg II A1239.925 is at 270 kms" 1 relative to N V A1238.810, 
and therefore can contaminate the LMC component of the strong 
line of the doublet (see panel 2 in Fig.[2j. 

No interstellar Galactic or LMC N V is detected towards any 
of the stars in our sample. The flux is much higher in the weaker 
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Figure 2. Same as Fig.[T]but towards Sk— 67°211. 

transition at 1242.804 A than in the transition at 1238.810 A, and 
therefore N V A1242.804 produces better upper limits. We report 
the 3a upper limit to the column density in Table [2] by assuming 
the line is on the linear part of the curve of growth. 



2.5.3 C IV and Si IV measurements 

For C IV and Si IV, the continuum placement and the measure- 
ments of v,N,b can be constrained by considering both lines of 
the doublet. The absorption profiles of these ions appear also to be 
separated into several components rather than the single broad ab- 
sorption observed in the O VI line profile. Furthermore the C IV 
and Si IV profiles reveal that in most cases the absorption of these 
ions is coupled, i.e. both ions show very similar profiles. These 
facts motivated us to use a Voigt profile fitting method. The fit re- 
sults and the errors were obtained using the Voigt component fitting 
software of Fitzpatrick & Spitzer (1997). The fit results listed are 
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Figure 3. Same as Fig. [T] but towards Sk— 69° 246 (the dotted line shows 
another possible choice for the continuum). 



based on a simultaneous fit to both members of the C IV doublet 
on one hand and the Si IV doublet on the other hand. However, we 
choose not to fit the C IV and the Si IV profiles simultaneously to 
avoid making an a priori assumption that these two ions trace the 
same gas. Rather, we prefer to show if this is the case (or not) a 
posteriori. For the STIS E140M instrumental spread function we 
adopted the profile from the STIS Instrument Handbook (Profitt et 
al. 2002). Data taken through the 0'.'2 x 0'.'2 slit have an line spread 
function with a narrow core and a broader halo. 

While we do not present the results of the AOD measurements 
for C IV and Si IV, we undertook them to cross-check the validity 
of the profile fitting results. We find very similar results between 
the AOD and profile fitting methods for the MW and IVC compo- 
nents. For the HVC component, the column densities of C IV and 
Si IV could be 0.05-0.2 dex higher with the AOD measurements 
because the velocity interval over which the profile was integrated 
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Table 4. Temperatures and column density ratios 



VI A1031 

Mg II A1239 (IMC) 



N V A1238 



N V A1242 





C IV A 1550 




200 
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Figure 4. Same as Fig.Q]but towards Sk-71°45. 



was larger than the breadth of the model. Hence, the profile fits to 
the HVC component are uncertain, and the AOD column density ra- 
tios (see i)2.6t are more reliable since there is no assumption on the 
profile of the absorption line. For the strong LMC components, the 
AOD provided in general a lower limit in agreement with the profile 
fitting measurements; and in the cases where the LMC components 
were essentially resolved (C IV and Si IV towards Sk— 65°22 and 
C IV towards Sk-69°246), the AOD and profile fitting gave the 
same results within the la errors. 

Below, we provide further insights on the choices of the con- 
tinua and component fitting for each line of sight. 

Sk—65°22: The C IV stellar wind is really strong, reduc- 
ing subsequently the flux near the MW component of the C IV 
A1548.195 (see Figs. [TJ and [5}. For Si IV, the chosen continuum is 
somewhat intermediate between plausible slightly higher and lower 
continua, but within the la errors the results can be reproduced 
by either solution. The MW, IVC, and HVC absorption can each 



Sight line 



T 

(10 5 K) 



N(C IV) /N (Si IV) 7V(0 VI) /N(C IV) 



LMC - H II region/Superbubble 
Sk-65°22 ^0.27 1.6 ± 0.3 

Sk-67°211 (<)0.2±0.4 7.1 ± f;§(:) 

Sk-69°246 (<)1.6±0.7 1.1±0.2(:) 
Sk-71°45 (<)0.2±0.3 2.4 ± 0.3(:) 







LMC - not associated with H II region 




Sk- 


-65°22 


< 4.2 


2.5 ± 0.3 


3.6 ±0.5 


Sk- 


-67°211 


6.0 ±3.0 


1.9 ±0.3 : 


2.5 ±0.3 : 


Sk- 


-69°246 


< 10 


1.5 ±0.5 : 


3.4 ± 1.4 : 








HVC 




Sk- 


-65° 22 


< 4.4 


3.3 ± 1.0 


7.4 ± 1.6 


Sk- 


-67°211 


< 12 




2.5 ±0.4 


Sk- 


-69°246 


1.5 : 


3.7 ± 0.9 


2.9 ± 1.7 


Sk- 


-71°45 




3.2 : ±0.7 : 


5.9 ±0.7 








IVC 




Sk- 


-65° 22 


< 2.2 


7.6 ± 0.6 


1.3 ±0.1 


Sk- 


-67°211 


< 3.8 




2.3 ± 1.6 


Sk- 


-69°246 


< 1.6 


3.9 ±0.6 


1.4 ±0.4 


Sk- 


-71°45 


4.7 ± 1.4 


4.4 ±0.4 


1.7 ±0.1 








MW 




Sk- 


-65° 22 


0.8 ± 0.4 


2.5 ± 0.3 


3.2 ±0.4 


Sk- 


-67°211 


< 1.9 




1.2 ±0.5 


Sk- 


-69°246 


1.5 ± 0.9 


3.5 ± 0.3 


1.6 ±0.4 


Sk- 


-71°45 


0.4 ±0.9 


2.6 ±0.4 


2.0 ±0.2 



Note: The temperatures are derived from the b-values of C IV and Si IV. 
The ";" symbol indicate some uncertainty in the result that is not reflected 
by the formal errors deduced from the measurements. The "(<)" symbol 
indicates that these lines are strong and may be composed of several 
components. The "(:)" symbol indicates that these lines are saturated and 
therefore the ratio is more uncertain than the formal errors from the profile 
fitting suggest, a: No O VI is found to associated with the narrow C IV 
component. See also £|2.6l for the high ratio measurements. 

be fitted well with one component. The HVC kinematical struc- 
ture remains largely uncertain and unknown; there is almost cer- 
tainly more than one component over the velocity range [100, 180] 
km 1 . The LMC absorption is fitted well with a broad and a nar- 
row component. 

Sk—67°211: Along this line of sight, the continuum to Si IV 
is uncertain, especially for the HVC component. Since there is a 
C IV HVC component, the continua to the Si IV lines are likely to 
be somewhat higher than seen in Fig.[2]in the HVC range, but there 
is no useful constraint to indicate how much higher. Our choice 
of continua therefore places a firm lower limit on the Si IV col- 
umn density of this HVC. A higher continuum for Si IV would also 
likely change slightly the IVC results. In contrast, the continuum 
to C IV is more straightforward. This is the only line of sight in 
our sample where Si IV and C IV in the MW and IVC components 
may trace different plasmas since the profile shapes and kinemat- 
ics of Si IV and C IV are not consistent with each other. For the 
LMC, the C IV can be fitted well with a broad and a narrow com- 
ponent. For LMC Si IV absorption, the fit is not unique: one can fit 
the LMC absorption with a narrow and a broad component or with 
three narrow components. We favor the two-component fit because 
a posteriori the kinematics match well the kinematics derived for 
C IV. This produces further uncertainties on the physical parame- 
ters derived for Si IV, and we highlight the uncertainty in Table [3] 
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by a colon (we discuss in Sj3]fhat a two-cloud model for the narrow 
C IV and Si IV components is likely). 

Sk-69°246: While for O VI this line of sight has the most 
complicated continuum placement near the interstellar O VI, for 
C IV and Si IV the continuum placement is comparatively straight- 
forward. For the MW, IVC, and HVC, we again use a single cloud 
model for each of these components. The C IV A1548.195 may 
suggest a two-component structure at about 3 and 25 kms -1 , and 
a narrower component for the IVC, but within the S/N and con- 
sidering the other transition, we did not manage to improve sta- 
tistically the fit with such a model. For the LMC, we used again 
two components. Yet a combination of two narrow or narrow and 
broad components can reproduce well the observed profiles. For 
C IV, the reduced-;^ 2 are essentially the same, 1.00 and 1.03 for 
the narrow-narrow fit and narrow-broad fit, respectively. For Si IV, 
the reduced-^ 2 is slightly better for the narrow-broad component 
fit compared to the narrow-narrow component fit, 1.12 compared 
to 1.27, respectively. The centroid of the broad component of Si IV 
and C IV matches well the centroid of the broad O VI absorption. 
The broadenings of O VI and C IV also appear compatible within 
la. Finally, the extra absorption observed in the red wing of C IV 
and Si IV can only be reproduced adequately with a broad compo- 
nent, putting more weight on this solution. We therefore favor the 
broad-narrow component fit, but we highlight the uncertainty in the 
broad component in Table[3]by a colon. We note that the LMC ab- 
sorption lines are very strong; if a narrower component is present 
in the profile, the column densities could be uncertain by an order 
of magnitude. 

Sk—71°45: This is the second giant star having a similar stel- 
lar spectrum near Si IV to that of Sk— 67°211. Although we had 
less difficulty finding continua that appear satisfactory for both 
lines of the Si IV doublet than for Sk— 67° 21 1, the true shape of the 
continua remains unknown and the quoted errors do not adequately 
reflect this systematic uncertainty. As for the other lines of sight, the 
continua for the C IV doublet lines are comparatively straightfor- 
ward to model. The C IV and Si IV absorption profiles are modeled 
with a single component for each cloud. For the HVC, we fixed 
the b-value for both ions; the centroid velocity for Si IV was also 
fixed to the C IV velocity. If b is not fixed, the fitting procedure 
tends to fit an unrealistically broad component. The LMC compo- 
nent of Si IV A1393.755 appears to have a wing not reproduced 
by the fit. However, this is not observed in the other transition, and 
we attribute this to the continuum placement. The addition of an- 
other LMC component to model this absorption fails to produce an 
adequate fit to the lines. The LMC lines are very strong and if a nar- 
rower component is present in the profile (and it is likely the case, 
see discussion in Sf3), the column densities could be again uncertain 
by an order of magnitude. 

2.6 Column density ratios of the highly ionised species 

Theoretical models predict the highly ionised species are produced 
via a range of mechanisms, such as in shocks or in interfaces be- 
tween different temperature gas, with conductive heating, turbulent 
mixing, or in radiative cooling gas (see references in Table [5] and 
summaries from Spitzer 1996, Fox et al. 2003 Indebetouw & Shull 
2004). O VI is of prime importance since in the galactic environ- 
ments, it is produced mainly via collisional processes. In contrast 
Si IV and C IV can be produced by photoionisation. If there is a 
kinematical relationship between the highly ionised species, the ra- 
tio of the high-ion column densities can be used to determine which 
physical mechanism is predominant in the ISM, if any, and how the 



mechanisms responsible for producing the highly ionised species 
may change in different galactic environments. 

Two methods were used to determine the ratios of highly 
ionised species. In Table [4] the ratios of the column densities of 
C IV and Si IV obtained from the profile fitting are summarised. 
We also measured systematically the O VI column densities and 
N V column density limits by using the velocity extent of each 
cloud defined by the profile fitting of C IV and Si IV. The ratios 
7V(0 Vl)/iV(C IV) derived from these measurements are listed in 
the last column of Table [4] Since the O VI absorption profiles do 
not break up into easily-identifiable clouds or blends of clouds like 
C IV and Si IV, O VI may or may not be associated with the other 
ions. 

In the second method, we compare the apparent column den- 
sity profiles over a certain velocity range that is listed in Table [5] 
This method does not make any physical assumption of the prop- 
erties of the gas. We separate the profiles in several velocity ranges 
corresponding to various clouds obtained in the profile fitting (MW, 
IVC, HVC, LMC). We give in Table|5]the range of ratios that spans 
the full range of the ±1<t errors. If only an upper limit is available, 
the 3cr upper limit and the apparent column density were estimated 
over the same velocity range. 

There is a good agreement between both methods for the MW 
and IVC components. This is not surprising since the MW and IVC 
components are strong, but not saturated, and well separated from 
one another and from other component blends. The component fits 
to the C IV and Si IV HVC profiles are not always unique due to the 
broad wing-like shape of the HVC absorption. Yet, there is a rough 
agreement between the AOD and profile fitting methods. Because 
the HVC profiles are not Gaussian, we favor the AOD results over 
this velocity range. Hereafter, we discuss separately the quite nar- 
row components of the highly ionised species in the LMC and the 
broad LMC components, which we assume includes all of the O VI 
given the lack of narrow absorption in this ion. 

For the LMC broad component, the ratios derived from the 
AOD and profile fitting methods can generally not be compared: 
the AOD estimates the ratios in the wing of the absorption where 
it is not contaminated by the narrow components observed in C IV 
and Si IV, while the profile fitting models the broad absorption even 
in regions where the strong narrow absorption is present (see Figs. [5] 
and|6]). In two cases, the comparison between the two methods can 
be made because the AOD estimates are made over a large frac- 
tion of the profile: Sk-67°211 («lsr = [180,230] kms -1 ) and 
Sk-69°246 (v LSR = [139, 190] kms^ 1 ). The AOD O VI/C IV 
ratios are in excellent agreement with those derived with the profile 
fitting for these sight lines. For the LMC narrow component, the 
only line of sight for which the ratios appear secure is Sk— 65°22; 
along this sight line the narrow component of C IV and Si IV is not 
too strong and, therefore, the profiles can be modeled without too 
much uncertainty. For the other lines of sight, the high-ion ratios 
of the narrow component are uncertain because C IV and Si IV are 
saturated. We discuss the high-ion ratios and implications in §4. 



3 OVERVIEW OF THE HIGH-ION KINEMATIC 
PROPERTIES 

3.1 Kinematics of the LMC absorption profiles 

Interstellar absorption from the low- and high-ionisation species is 
always detected over the full range of velocities from the Milky 
Way to the LMC (except for N V) along the four lines of sight in 
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Table 5. Ratio of the highly ionised species from AOD comparison 



Model/Sight line 


N(0 VI) /N (Si IV) 


N(0 VI) /N(C IV) 


N(0 VI)//V(N V) 




Models 






Collisional Ionisation Eq. (CIE) 


^ 111 


3.7-200 


0.65-55.6 


Radiative Cooling (RC) 


> 76.9 


5.9-27.8 


11.1-18.9 


Conductive Interface (CI) 


^ 17.2 


1.1-23.8 


2.0-15.4 


Turbulent Mixing Layers (TMLs) 


0.8-11.5 


0.06-1.0 


1.7-5.9 


Shock Ionisation (SI) 


^ 5.6 


0.96-62.5 


19.6-32.3 


Supernova Remnant (SNR) 


74.4 


5.2-8.4 


11.9-16.0 


Stellar Wind (SW) 


100 


6.3 


15.8 


Halo SNR 




0.42-5.0 


3.0-11.1 




LMC - H II rej 


;ion 




Sk-65°22 


< i-o 


£0.8 




Sk-67°211 


< 0.3 


< 0.1 




Sk-69°246 


< 0.8 


< 0.4 




Sk-71°45 


< 0.5 


<0.2 




LMC - not associated with H II region 




Sk-65°22 [187,232] 


> 37 


> 5.4 


> 3 


Sk-67°211 [180,230] 




2.2-5.1 


> 7 


Sk-69°246 [139,190] 


4.6-10.7 


1.7-5.0 


> 10 


Sk-71°45 [162,190] 




5.4-9.1 


> 9 


Sk-71°45 [260,320] 




> 11.0 


> 6 




HVC 






Sk-65°22 [88,187] 


20.4-57.5 


4.7-9.5 


> 5 


Sk-67°211 [100,182] 




2.94.3 


> 6 


Sk~69°246 [85,139] 


8.1-30.1 


1.4-5.6 


> 2 


Sk-71°45 [104,183] 




4.9-7.2 


> 9 




IVC 






Sk-65°22 [23,88] 


4.0-12.6 


0.4-1.3 


> 6 


Sk-67°211 [30,100] 


7.9-20.1 


0.9-2.2 


> 11 


Sk-69°246 [41,85] 


3.9-7.9 


0.9-1.7 


> 2 


Sk-71°45 [26,104] 


3.2-20.0 


1.0-3.2 


> 10 




MW 






Sk-65°22 [-45,23] 


2.5-11.2 


1.1-3.2 


> 4 


Sk-67°211 [-43,30] 


4.0-7.1 


1.0-1.8 


> 6 


Sk-69°246 [-61,41] 


3.2-6.3 


1.1-2.0 


> 2 


Sk-71°45 [-51,26] 


2.5-7.1 


1.0-1.8 


> 4 


Note: See ^12-61 for the high ratio measurements. References for the models: CIE 


- Sutherland & Dopita 


(1993) for gas temperature T = 


= (2-5) X 10 5 


K; RC - Edgar & Chevalier (1986); CI: 


iBorkowski, Balbus, & Fristroml Jl990l 


) for magnetic field 


orientations 9 = 0- 


-85° and interface gas 


aaes 10 5-7 vr: TML -ISlavin. Shull. 


& Beeelmanl J 19931) 


for a gas flow in the range 25-100 km s 1 


and temperature of postmixed gas in 


the range log T = 


5-5.5; SI -iDooita & Sutherland Jl996l) for 



j 19921) for SNR ages 10 5 
10 6.0-7.2 yr 



yr: 



SW - IWeaver et al] j 19771) : Halo SNR 4Sheltorj 09981) for SNR ages 



Table 6. Summary of the velocity centroids in the LMC 



Name 


v(Ha) 
(kms -1 ) 


v(Fe III) 
(kms -1 ) 


v{0 VI) 

(kms -1 ) 


v(C IV) (narrow) 

(kms -1 ) 


u(C IV) (broad) 
(kms -1 ) 


Environments" 


Sk-65°22 
Sk-67°211 
Sk-69°246 
Sk-71°45 


270 ± 15 
285 ± 17 

252 
227 ± 19 


270.6 ± 1.3 
286.3 ± 1.0 

214.7 ±3.1" 

214.8 ± 1.8 


258.9 ± 1.7 

251.2 ±2.3 
215.7 ±2.1 

227.3 ± 1.2 


265.4 ±0.5 
272.2 ±0.3 
214.1 ±0.5 
221.0 ±0.3 


271.7 ±1.7 
251.2 ±3.4 
213.5 ± 1.8 


Faint H II region 
H II region 
Superbubble (X) 
Superbubble 



Note: a: For this sight line, Fe III was fitted with a two-Gaussian component that yields for the second component 
«LSR = 242.3 ± 4.4 kms -1 . 
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Figure 5. Normalised profiles of the selected ions observed along each line of sight. In the S II panel, the gray histogram represents the O I A 1302 absorption 
profile. The solid black curve shows the profile fit to the C IV and Si IV lines. We note that the continuum placement near Fe III is typically uncertain, and 
these profiles are shown to demonstrate only the velocity distribution of the Fe III absorption. The vertical dashed line represents the systemic Ha velocity. 



our sample. In Figs.[5]and[6] we show the normalised profiles of the 
various species. In these figures, the vertical dashed line represents 
the profile-weighted average velocities measured for the Ha pro- 
files (results from Danforth 2006 and Chu et al. 1994, see Table[6]l. 
Both superbubbles show evidence for expanding shells in Ha; the 
projected radial velocities of the approaching side are shown by the 
dotted lines in Fig.[6](Y.-H. Chu 2006, priv. comm.; see also Chu & 
Kennicutt 1994). There is no obvious structure of expanding shell 
seen in Ha toward the other sight lines. The average velocities ob- 
tained from the Ha emission profiles agree remarkably well with 
the centroids measured in the Fe III profiles (see Table|6]l. The Fe III 
profiles in the spectra of Sk— 69°246 and Sk— 71°45 show indica- 
tion of two components with one matching the systemic velocity of 



the Ha profile and the other matching the velocity of the expansion 
of the shell. The agreement between the Ha and Fe III profiles is 
quite remarkable given that Ha emission data probe gas in front of 
and behind the star while absorption spectra only probe gas in front 
the star. 

Figs. [5] and [6] show continuous absorption at LMC veloci- 
ties (~130-310 kms -1 towards Sk-65°22; ~210-370 kms" 1 
towards Sk-67°211; —130-310 kms" 1 towards Sk-69°246; 
~150-340 kms" 1 towards Sk-71°45) for O VI, Si II, and O I. 
(We note that velocities 80 kms" 1 < ulsr < 130 kms -1 may 
have an LMC origin as well; see §4.2.) Over these velocities, 
the O VI profiles are broad and have vaguely Gaussian shapes. 
The weakly ionised species consist of several narrow components. 
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Figure 6. Same as Fig.[3]for the two sight lines Sk— 69°246 and Sk— 71°45. Towards Sk— 69°246, we show the normalised profiles of O VI from both the 
high and low continua). The S II panels also include the O I A 1302 profiles. The vertical dotted line represents the expanding Ha shell velocity (Y.-H. Chu 
2006, priv. comm.). 



Since O VI traces much hotter gas than Si II or O I, the difference 
in the appearance of the profiles is expected. In contrast with the 
weakly ionised species, the overall kinematic morphology of the 
O VI profiles is quite similar between the various sight lines. 

Si II and O I also reveal the presence of an IVC with respect to 
the LMC at i>lsr ~ 150-220 kms -1 , clearly separated from the 
main absorption by roughly —50 kms -1 . This LMC-IVC com- 
ponent possibly traces some ejecta from the LMC disk (see §4.2). 
The O VI absorption is always present in this IVC component. We 
note here that the important difference between the choice of strong 
transitions compared to the weaker Fe II absorption profiles used by 



lHowketal.N2002h : the strong transitions reveal cooler, neutral and 
partially ionised gas traced by weakly ionised species and hot gas 
traced by O VI are systematically found at similar velocities. 

The morphology of the C IV and Si IV profiles differs from 
the morphology of the O VI profiles. Along our four sight lines, 
a strong narrow component is observed in C IV and Si IV but not 
revealed in O VI. Along three sight lines, the profiles of C IV and 
Si IV have the strong narrow component blended with a broad com- 
ponent (see i)2.5.3l >. The broad component is likely associated with 
some of the observed O VI. There is a good agreement between the 
centroids of broad C IV and Si IV components with the centroid 
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of the broad O VI absorption towards Sk-67°21 1 and Sk-69°246 
(see summary in Table [6}. Towards Sk— 65°22 the broad compo- 
nent of C IV and Si IV is shifted by 12 kms -1 because the lower 
velocities observed in the O VI profile are not seen in the C IV and 
Si IV profiles. 

In the faint H II region where Sk— 65°22 is located, the cen- 
troids of the narrow components of C IV and Si IV align with those 
of Fe III and Ha as well as with the main component of the weaker 
ions (see Fig. [5}. This sight line has the simplest LMC kinematics 
of our four. For Sk— 67°211, the narrow components of the C IV 
and Si IV absorption profiles appear shifted to lower velocities rel- 
ative to the velocities of Fe III and Ha. However, C IV and Si IV 
are saturated and the narrow profiles may contain more than one 
component. Towards Sk— 69°246, the main absorption of C IV and 
Si IV is observed at the expanding velocity seen in Ha and Fe III. 
There is very little absorption in the C IV and Si IV profiles at 
252 kms -1 , which corresponds to the systemic velocity of Ha. Fi- 
nally, towards Sk— 71°45, also located in a superbubble, the C IV 
and Si IV profiles do not align with either component of Fe III, 
their velocities being somewhat intermediate between the systemic 
and expanding shell velocities observed in Fe III and Ha. The gas 
traced by the narrow components of C IV and Si IV, and Fe III 
are kinematically related in a different manner toward each sight 
line: there is an excellent agreement in their velocity profiles in the 
quiet H II region; in the bright H II region, C IV and Si IV is blue- 
shifted with respect to Fe III, possibly revealing an highly-ionised 
expanding shell; in the X-ray bright superbubble, C IV and Si IV 
appear to probe mostly the expanding shell seen in the Ha profile; 
and in the superbubble DEM 221, C IV and Si IV do not appear to 
directly be kinematically related to Ha, with the velocities of the 
highly ionised species intermediate between the expanding or sys- 
temic velocities of Ha. We finally note that weaker ions (Si II and 
O I for example) and higher ions (O VI) can produce absorption at 
these velocities. 

In summary, the absorption profiles of the weakly ionised 
species and O VI always span the same velocity range, though 
the former are always found in discrete, narrow components, while 
the latter is broad and widespread. There is either a positive, nega- 
tive, or no shift between the average velocities of O VI and weakly 
ionised species. However, the comparison of the average velocities 
of O VI with the weakly ionised species is complicated since O VI 
is almost certainly a blend of several components (see §3.2). The 
velocity extent of the C IV and Si IV LMC profiles is sometimes 
smaller than the O VI profile extent. The narrow LMC components 
of C IV and Si IV may be kinematically connected to the gas traced 
by Ha and the weakly ionised species (e.g., Fe III, Si II), but not 
in a systematic way, implying that even though the narrow compo- 
nents of C IV and Si IV probe cool gas (T < 2 x 10 4 K, see below) 
as Fe III, their origin/production may differ, at least in some cases. 

3.2 Broadening of the high-ion absorption profiles 

The gas temperatures derived from the 6- values of C IV and Si IV 
are given in Table [3] These remain quite uncertain for the LMC 
narrow component as the profiles are very strong, and the errors 
do not reflect the possibility that the profiles have more than one 
component. To derive these temperatures, we assume that C IV and 
Si IV reside in the same gas, which is supported by the compati- 
bility in the centroids of C IV and Si IV for most clouds. In Ta- 
ble [3] we list the temperature of the gas derived from the observed 
broadening of the absorption lines. If b(C IV) > 6(Si IV) then T 
can be determined from b = y/2kT/(Amn) + b 2 nt , where Ont IS 



the non-thermal component and the other symbols have their usual 
meaning. If b(C IV) < 6(Si IV) an upper limit on the temperature 
can only be derived because we cannot estimate the non-thermal 
broadening. 

For the strong, narrow LMC component (denoted LMC (H II) 
in Table [3}, we find a good agreement between the centroids of 
Si IV and C IV. The Doppler broadenings suggest a temperature of 
2 x 10 4 K along three lines of sight (towards Sk-69°246, within 
2<j error, the temperature is compatible with 2-9 x 10 4 K). Since 
the profiles are strong and saturated towards at least Sk— 67°211 
and Sk— 69°246, the derived b is uncertain but places a firm upper 
limit on the temperatures; these clouds likely have T < 2 x 10 4 
K. The non-thermal broadenings of C IV and Si IV appear to be 
as large as the thermal broadening if the profiles consist of a single 
component. 

There must be little or no O VI associated with the narrow 
and strong C IV and Si IV components since there is no indication 
of a narrow component in the overall broad and smooth profile of 
O VI (see the limits of O VI to C IV or Si IV ratios in Table |5). 
Hence there is no evidence for O VI arising directly within the 
superbubble shells or the H II regions. We also note that these H II 
regions and superbubbles do not produce any detectable N V. This 
is not surprising since the temperatures deduced from the C IV and 
Si IV line-width are only a few x 10 4 K or less, too cool to produce 
any significant amount of O VI or N V0 

For the broad C IV and Si IV components of the LMC, their 
6- values imply temperatures of a few x 10 s K or less. Their broad- 
ening is consistent with that of O VI, and these components likely 
probe hot gas. For O VI only upper limits to the temperatures can 
be derived, and these are always consistent with temperatures of 
a few 10 5 K. The broad LMC components of the highly ionised 
species have fe(LMC) ^ i/6 2 (MW) + b 2 (IVC), suggesting that 
they likely trace more than two different clouds. 



4 THE HIGHLY IONISED PLASMA IN THE LMC 

We found that the absorption of the highly ionised species at veloc- 
ities similar to the LMC disk is dominated by a narrow absorption 
observed in C IV and Si IV but not in the O VI. The O VI along 
all four sight lines is broadly distributed, and the profile fitting of 
C IV and Si IV reveals broad absorption in these ions along 3 of 4 
sight lines. These two types of absorption components are certainly 
produced by different mechanisms. We discuss them separately. 

4.1 The narrow absorption components 

The photon flux from hot stars drops sharply above 54 eV because 
of the He II ionisation edge, so photoionisation by hot stars can 
produce Si IV (33.5 eV) and C IV (47.9 eV) but not significant 
N V (77.5 eV) or O VI (1 13.9 eV). This e nergetic argument was, 
for example, advanced by iFox et al.l |2003) to explain the narrow 
C IV and Si IV absorption components observed towards a Galac- 
tic sight line. We show in the previous section that the broadenings 

1 N could be underabundant in the LMC H II regions (Garnet 1999, Welty 
et al. 1999). However, the ratios of C IV/N V for the narrow compo- 
nent are greater than 1.7, 478, 7.6, 10.0 towards Sk-65°22, Sk-67°211, 
Sk— 69° 249, and Sk— 71°45, respectively. So even if N may be intrinsically 
underabundant by a factor up to 6, with such high ratio limits (especially to- 
wards Sk— 67° 21 1), the bright H II region and superbubbles do not produce 
any observable amount of N V. 
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of the C IV and Si IV imply temperatures of a few 10 4 K and likely 
less in the stronger features because of unresolved components. The 
velocities of the narrow C IV and Si IV components are consistent 
in two cases with the velocities of the H II region or the expand- 
ing shell as traced by Ha emission and Fe III absorption. These 
observed properties are in agreement with photoionisation playing 
an important role in the production of these ions along with Fe III 
and Ha. But in the two other cases, the bulk of the absorption in 
C IV and Si IV do not correspond to either the systemic or expand- 
ing velocities of the regions probed, suggesting that additional or 
different production mechanisms are in play. 

The C IV/Si IV ratios in the narrow components appear al- 
ways greater than 1 (from AOD and profile fitting comparisons), 
although we remind the reader that the Si IV continua in the spec- 
tra of Sk-67°211 and Sk-71°45 are uncertain (see 323731 . For 
Sk-67°211, Sk-69°246, and Sk-71°45, the hidden complexity 
within the strongly-saturated C IV and Si IV profiles are likely to 
have more than one component, and our column densities are thus 
uncertain. While the O VI/C IV and O Vl/Si IV ratios are uncertain, 
they appear to always be smaller than unity, implying very little 
O VI is produced by the same mechanisms that produce the strong, 
narrow components of C IV and Si IV. This is further strengthened 
by the high C IV/N V ratios (see footnote[T). 

Photoionisation models of diffuse H II regions generally pro- 
duce ratios C IV/Si IV le ss than unity, although the addition of X- 
rays can boost this ratio l Black et al Ill98d :lc owie, Taylor. & Yorkl 



1 198 it iKnauth et alj|2003l) . The ft-values of the narrow C IV and 
Si IV components seen in our data strongly suggest, however, that 
photoionisation may play a role in producing this material. We have 
undertaken simple photoionis ation models using t he Cloudy code 
(v06.02.09; last described by iFerland et all 1 19981) to test the as- 
sumption that the underlying structures may give rise to the nar- 
row components. We calculated a model spherical H II region 
about a single star having properties simila r to the star Sk— 6 5°22: 
log L/L Q = 6.2, T cff « 40000 K (fr omlMassa et alj|2003h. W e 
adopted LMC metal abundances from iRussell & Dopital 1 1992b . 
Our models include a treatment of interstellar grains for heat- 
ing/cooling and opacity purposes; this dust affects the abundances 
due to elemental depletion (see Cloudy documentation). We as- 
sumed a clearing within the inner-most parsec of the H II region and 
constant density outward of this radius. In general we find that such 
models, with an extremely bright, high temperature star, can repro- 
duce (or exceed) the measured C IV/Si IV ratios observed in the 
narrow components along the sight lines in this work. We examined 
models without the depletion (i.e., a non-self consistent model with 
grains for thermal and opacity effects, but leaving all the metals in 
the gas phase), which verified the depletion effects did not domi- 
nate the C IV/Si IV results. These single star H II region models, 
while able to match or exceed our observed C IV/Si IV measure- 
ments, do not simultaneously match this ratio and the total column 
densities of the highly ionised species seen towards Sk— 65°22. In 
fact, the models that produce a sufficiently high column density 
tend to over-predict the C IV/Si IV ratio. However, these are ex- 
tremely simple models that do not take into account the other stars 
in the vicinity of our target objects, the complex geometries and 
density structures within the regions, etc. Given the high temper- 
atures and the quite high luminosities of our target stars, as well 
as their locations within star forming regions in the LMC, it seems 
likely that direct photoionisation from OB-type star photospheres 
could explain the narrow C IV and Si IV along some sight lines. 

Given the kinematic complexity of the other sight lines (see 
Sj3}, several processes may be at work or that the structures giving 



rise to the absorption have complex internal kinematics. In partic- 
ular, early-type st ars with strong stellar wind create large cavity in 
the ambient ISM dCastor. McCrav. & Weavedl 19751) . Bubbles have 
complex kinematic s and have a mixtur e of hot (T > 10 6 K) and 
cool (T < 10 4 K) JWeaveretai]|l977l) . The hot gas is principally 
found in the shocked wind. The cooler gas mainly is found in the 
shell behind the conducting front following a rapid radiative cool - 
ing (see also. lShapiro & Moord[~976l ; ISufherland & Dopitdll993l) . 
Such process also creates an interface between the hot and cool 
gas at T ~ 10 5 -10 6 K where some of the observed O VI is pro- 
duced (see below). The ionisation of t he outer shel l depends mostly 
on the radiation of the central star l lWeaver et alj|l977h . It is in- 
teresting to note that the two hotter stars in our sample show the 
strongest C IV and Si IV absorption associated with an expanding 
shell (see Sj3)- The stellar wind model therefore provides a qualita- 
tive explanation for part of observed narrow component associated 
with the outer/expanding shell. This is also supported by the re- 
cent Chandra observations of 30 Do, which shows that Sk— 69° 246 
may be a collidi ng-wind binary with bright X-ray emission from 
2.1 keV plasma dTownslev et alj2 006). However, 30 Dor and other 
regions of the LMC are also well known for their superbubbles cre- 
ated not only by stel lar winds but also SNRs. The SNR models of 
ISlavin & Coxl (l992) can also produce cool highly ionised at 10 
K. This is not surprising since the SNR and stellar wind models ex- 
plore very similar physics; in the former the source of the shocked 
is the supernova, in the latter case, it is the central star. The stellar 
wind and SNR models place the expanding shell within a radius 
of about 30 to 100 pc from the central region. Yet, other physical 
scenarios might produce such cool highly-ionised gas, as well. For 
example, emission of a hot, cooling plasma could photoionise the 
gas; O VI and N V are negligible with respect to C IV and the 
temperature of the gas is low (T < 10 4 K ) so the C IV and Si IV 
broadenings are small (Knauth et al. 20q3). 

In summary, in view of the kinematical differences between 
the various sightlines, the origins of the narrow components ob- 
served in the C IV and Si IV profiles are most likely multiple. In 
the spectra of Sk-67°211 and Sk-71°75, it is likely that the nar- 
row components of C IV and Si IV are in fact at least two narrow 
components (hence tracing gas at T < 10 4 K) at the systemic and 
expanding shell velocities. In the spectrum of Sk— 65°22, the ab- 
sorption is comparatively weak and consistent with a single nar- 
row component associated with the H II region around the star. 
And in the spectrum of Sk— 69°246, there is only evidence of a 
strong absorption at the velocity of the expanding shell, showing 
the peculiarity of this line of sight. While the exact origins of the 
highly ionized species probed by the narrow components remain 
unknown, their properties can be chiefly associated with the stellar 
environments. 



4.2 The broad absorption components 

Towards each line of sight, the O VI LMC component is very broad 
with b ~ 50 km s ~~ 1 , much larger than the broadening of the Galac- 
tic plus the IVC or the temperature at which O VI peaks in abun- 
dance in CIE, implying that the LMC O VI component is composed 
of several kinematic components. There is evidence for broad C IV 
and Si IV LMC absorption components associated with the O VI 
towards at least 3 lines of sight. Although the velocity structure 
of the broad absorption is not fully known, the high-ion profiles 
suggest that this gas traces hot collisionally ionised gas at temper- 
atures of a few times 10 5 K. We also showed in Sj3]fhat O VI and 
Si II have absorption over nearly the entire range of velocities from 
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Figure 7. Summary of the column density ratios of VI to C IV for the 
MW, IVC, HVC and broad LMC components towards the four sight lines, 
obtained from the AOD measurements (see Table|3J. 



—50 to +330 kms" 1 . C IV and Si IV cover a smaller velocity 
range. In particular they do not show absorption at the lowest or 
highest velocities associated with the LMC. Thus, gas at the lowest 
and highest velocities makes a proportionately smaller contribution 
to the total columns of C IV and Si IV than it does in O VI. 



4.2.1 Ratios of the highly ionised species 

Fig. [7] shows the O VI/C IV ratio (from Table [5} along our four 
sight lines, separated into the four principal components: the MW 
(wlsr < +30 kms" 1 ), IVC (30 < v LSR < 100 kms" 1 ), HVC 
(100 < -ulsr < 175 kms" 1 ), and broad LMC (w LS r > 175 
kms" 1 ) velocity components. Over the velocities where the nar- 
row components discussed above are present the ratios are in- 
ferred mostly from the profile fitting (Table [4}, except towards 
Sk— 67°21 1 and Sk— 69°246 where the absorption is weak enough 
that the ratios of the apparent column densities can be estimated 
over a portion of the absorption profiles. In the broader wings of 
the absorption profiles, the ratios are inferred from our AOD anal- 
ysis (TableO- For Sk— 65°22, the fits to the broad components of 
C IV and Si IV absorption appear secure, but the fits to the broad 
absorption along the other sight lines are not unique (see i]2.5.3l >. 
We found a good agreement between the ratios estimated from the 
AOD and profile fitting for Sk-67°21 1 and Sk-69°246, strength- 
ening the results from the profile fitting. 

While there is some room for variations in the O VI/C IV ratio 
as a function of velocity along a given sight line (our approach to 
deriving the values plotted in Fig.|7]is relatively conservative), the 
HVC and broad LMC components have quite similar O VI/C IV 
ratios, and the MW and IVC components are generally in agree- 
ment with one another, but smaller than the HVC-LMC values. 
The observed O Vl/Si IV ratios for the LMC and HVC are roughly 
consistent with those found in the MW or IVC components, but 
the O VI/C IV ratios are generally quite a bit larger in the HVC 
and LMC material than in the MW or IVC material. Moreover, the 
O VI/C IV ratios in HVC and LMC gas seem to vary between sight 
lines, but the ratios in the HVC and LMC components along a given 
sight line are always very similar. 



4.2.2 Evidence for outflows 

The similarity of the O VI/C IV ratios between the HVC and LMC 
gas suggests that (i) the HVC component has its origins in the 
LMC, and (ii) the processes producing the highly-ionised gas in the 
LMC affect a large range of velocities along a given sight line in 
a similar man ner. The first of these conclu sions supports a conclu- 
sion drawn bv lStavelev-Smith et alj J2003h in their study of the H I 
structure of the LMC. These authors found that the relatively high 
column density clouds at velocities ulsr ~ 100-160 kms -1 to- 
wards the LMC are often seen projected onto H I voids in the LMC 
disk, and they are connected to the disk with spatial and kinematic 
bridges in position-velocity plots. They also note the presence of 
material at velocities outside the range of velocities of the LMC 
disk. They conclude based upon these lines of evidence that the gas 
in these clouds is likely expelled from the LMC disk. 

Our absorption line measurements are much more sensitive 
to low colu mn density matter than the H I emission line mea- 
surements of lStavelev-Smith et al. I J2003h . Thus, we are potentially 
probing gas with significantly different ionisation states. Indeed, 
Staveley-Smith et al. only display maps for the HVC material with 
N(H I) > 10 18 cm -2 , which shows most of the gas concentrated 
to the southwestern edge of the LMC. Our lines of sight are con- 
centrated significantly further to the east than most of their clouds. 
iDanforth et alj d2002h , however, have demonstrated that HVC ma- 
terial can be found over a substantial fraction of the LMC using 
FUSE absorption line measurements. Their measurements of Fe II 
equivalent widths towards LMC stars show a broad swath of high- 
velocity material covering the eastern edge of the galaxy, with the 
highest equivalent widths arising in the northeast quadrant. (Note 
that the coverage of the southwestern quadrant in their data was 
minimal.) 

The most promising way to confirm our conclusion that the 
HVC material is associated with the LMC is to obtain absolute 
abundances of the HVC gas, which will require more sensitive H I 
emission line observations towards these sight lines. An absolute 
abundance derived from the ratio of O I to H I similar to the abun- 
dance of the LMC would strongly suggest an LMC origin. 

While we acknowledge the uncertainties in this conclusion, 
we nonetheless proceed assuming that the HVC towards the LMC 
is associated with this galaxy. Thus, the velocity range for LMC 
gas along these sight lines should be defined as everything at 
«lsr > 100 kms -1 . The velocities of the LMC IVC and HVC 
components indicate that they are expanding away from the warm 
disk gas or, for the gas at lower absolute velocities relative to the 
LMC disk, that the gas is extended above the disk, perhaps rotat- 
ing in a thick disk. Only towards Sk— 71°45 is there evidence for 
both significantly red-shifted and blue-shifted components in the 
weakly ionized species and O VI. Thus, little of the material we see 
can represent gas returning to the LMC disk. The velocity of the 
HVC relative to the LMC disk is larger than a simple estimate of 
the escape velocity of the LMC, v e sc ~ V2v c j rc ~ 1 00 kms -1 
(with v C i rc — 72+7 kms" 1 from lAlves & Nelsonl2000l) , implying 
that the material may escape the LMC, polluting the intergalactic 
space between the LMC and the Milky Way or serving as fuel for 
the Magellanic Stream. 

The LMC-IVC may trace an interstellar thick disk about the 
LMC, perhaps similar to tha t seen about the Milky Way (e.g., 
ISavage, Sembach. & Lul 19971) . However, it seems unlikely that the 
HVC material can be associated with an extended, dynamically- 
supported thick disk component due to its extremely large veloc- 
ity separation from the LMC. Theoretical modeling and stellar 
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observations of the LMC have revealed the existence of a stellar 
thick disk (e.g., I Weinberg! I2OO0I ; iBekki & Chiball2005l ; ICole et all 
2005). Neither the simulations nor the observations of evolved stars 
demonstrate velocity dispersion greater than about 20-40 kms" 1 . 
While thick disk gas may not follow the stellar kinematics, it is un- 
likely that the velocity dispersion of the gas would be larger than 
that of stars since the gas is more subject to hydrodynamical effects 
than stars (i.e. gas is not collisionless). Gas corotating in a thick- 
ened disk may appear at lower velocities due to projection effects 
asso ciated with the inc lination of the LMC disk to our line of sight 
(see lHowk et~al]|2002l) . While this could contribute to the velocity 
shift of the LMC-IVC relative to the thin disk material, it does not 
produce a large enough velocity shift to explain the HVC veloci- 
ties. Thus, the HVC material seems unlikely to be associated with 
a thick disk of material about the LMC. It is more likely tracing a 
large-scale outflow of material. 

4.2.3 Physical origins 

The comparison of the column densities of the highly ionised 
species can also be used to unravel their origins. The mode ls sum- 
mariz ed in Table [5] can be divided in 3 main mechanisms (Spitzer 
1996): (i) conductive heating mechanisms (including the SNR, stel- 
lar winds, conductive interface, and shock ionisation) in which cool 
gas evaporates into adjacent hot medium (the highly ionised gas at 
T ~ a few x 10 being found at the interface between the cooler 
and hotter plasmas); (ii) radiative cooling models in which hot gas 
cools, (iii) the turbulent mixing layers where hot gas mixes with 
cool gas. The broad component of the highly ionised species there- 
fore always probes indirectly the presence of hot gas (T > 10 6 K). 
In real physical environments, it is likely that these mechanisms are 
actually mixed all together to a certain degree, but depending on the 
physical conditions, one may dominate the other. 

We discussed in the previous section that strong stellar winds 
(or SNRs) most certainly occur, at least in the bright H II and su- 
perbubble regions. Therefore some of the O VI and broad C IV 
and Si IV absorption must originate in the interface between the 
shocked plasma and the swept-up interstellar gas. In velocity space, 
most of this absorption should occur near the expanding shell ve- 
locity where part of the strong narrow components C IV and Si IV 
are observed. While the observed high-ion ratios are very roughly 
consistent with those models, the observations also suggest that the 
environments very local to our target stars have little discernible 
effect on the O VI/C IV ratios of the broad absorbing components. 
The O VI/C IV ratios are the highest towards the faint H II re- 
gion associated with the supergiant shell DEM 48 (Sk— 65°22) and 
towards the superbubble DEM 221 (Sk-71°45). They are the low- 
est towards the H II region DEM 241 (Sk-67°211) and the X-ray 
bright superbubble 30 Dor (Sk— 69° 246). This difference between 
the two superbubble sight lines and the two H II region sight lines 
makes it impossible to conclude that one type of structure is respon- 
sible for producing higher or lower O VI/C IV ratios. The absence 
of clear signature in the high-ion ratios may be because the indi- 
vidual clouds cannot be distinguished. A very large fraction of the 
observed O VI (and broad component of C IV and Si IV) must also 
have a different origin than the stellar environment itself (stellar 
wind and SNR models predict N(0 VI) ~ a few x 10 13 cm" 2 , 
much less than the observed total O VI column density) and is not 
spatially localized near these H II regions and superbubbles even 
though some of it overlaps these structures in velocity space. 

When the observed LMC-HVC high-ion ratios are compared 
to theoretical models (see Table O, the turbulent mixing layer 



(TML) model always fails to produce large enough O VI/C IV; this 
is at odds with gas at the Milky Way velocities for which TMLs 
seem to be consistent with the observed ratios of highly ionised 
species (see Tables [4] and O. One caveat to this conclusion is that 
if the oxygen abundance is systematically enhanced with respect 
to carbon, it would be likely that more O VI is observed than in 
"typical" regions in the Milky Way. But enhancement of oxygen or 
deficiency of carbon are not observed in the LMC (see Appendix 
in Welty et al. 1999, and references therein). 

Another systematic between the various sightlines is deduced 
from Figs. [5] and |6] (see also discussion in Sj3j: cooler gas traced 
by Si II and O I is observed at all velocities where O VI absorp- 
tion is detected. The close correspondence in the total velocity ex- 
tent of the observed Si II and O VI absorption strongly suggests 
that the high and the weakly ionised species are kinematically cou- 
pled. A correspondence between O VI and low-ion velocity extent 
has a l so been noted in th e Milky Way ( Ho wk. Savage. & Se mbach 
2003;; iBowen et alj|2006T) . Such a relationship is expected if much 
or all of the O VI arises in interfaces between warm and hot gas. 
iBorkowski. Balbus. & Fristroml (l990) present models of conduc- 
tive interfaces in such an arrangement; they predict that if the O VI 
is found in both the evaporative and conductive stage, the O VI 
column density would be less than a few 10 13 cm -2 , implying 
that each line of sight in our work intercepts 5-10 interfaces. The 
broadening of the high-ions is consistent with such an interpreta- 
tion since, if there are 5-10 interfaces, each component of O VI 
could still have a temperature of a few times 10 5 K. While the low 
ion profiles are mostly unresolved, their velocity extents are consis- 
tent with a large number of components, and therefore with large 
number of interfaces. The profiles of weakly and highly ionised 
species cover a large enough velocity range that they are likely not 
simply tracing components associated directly with the H II regions 
or superbubbles housing the target stars. Thus, multiple interfaces 
are certainly possible and likely associated with gas out of the disk 
of the LMC. If conductive interfaces are the main origin for the 
broad component of the high ions, they have to be at least 3 x 10 
year old to reproduce the high-ion ratios. We note that the observed 
O Vl/Si IV ratios are somewhat smalle r than the ratio predicted by 
IBorko wski. Balb us. & Fristroml dl990h . but Si IV can be easily pho- 
toionised by other sources (e.g., photons from stars or the absence 
of the incorporation of He II A304 photons in the interface model 
that are capable of ionising efficiently Si III see Savage, Sembach, 
& Cardelli 1994), and this is not taken into account in their models. 
The limits derived from N V are always consistent with the con- 
ductive interface models. It is not clear within the frame work of 
the conductive interface models why the high-ion ratios would be 
similar over the large velocity range between the LMC and HVC 
components, and yet so dissimilar between the observed sightlines. 

In view of the complexity of the profiles and of the physical 
regions probed, other models may be in play, especially in the event 
that the apparent kinematic coupling does not imply coincidence of 
the highly and weakly ionised gas: 

- Strong shocks with u s hock > 200 kms" 1 would produce 
enough O VI to explain the observed ratios. Such strong shocks 
have not been inferred in the denser pa rts of H II and superbub- 
ble regions of the LMC dOev etal.l 2000). But high- velocity shocks 
are likely to arise in low-density gas where they remain undected 
since they would not contribute to the nebular emission (Y.-H. Chu 
2006, priv. comm.). The strong shock models are attractive because 
they provide a straightforward method for producing a consistent 
O VI/C IV ratio over a large range of velocities. If such strong 
shocks occur in the LMC, the gas traced by the LMC-IVC/HVC 
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Table 7. Interstellar O VI and environments in the LMC 



Id. Star R.A. Dec. logJV(OVI) Environments 

(°) (°) (dex) 



1 


Sk-67°05 


72.58 


-67.66 


13.89 ±0.07 


H II region DEM 7 


2 


Sk-67°20 


73.88 


-67.50 


14.26 ±0.10 


Field star, no Ha 


3 


Sk-65°22 


75.34 


-65.87 


14.14 ±0.02 


Supergiant shell DEM 48 


4 


Sk-66°51 


75.79 


-66.68 


14.31 ±0.07 


H II region DEM 56, no Ha 


5 


Sk— 67 69 


78.58 


-67.13 


14.48 ± 0.03 


H II region DEM 107 


6 


Sk-68°80 


81.62 


-68.84 


14.61 ±0.04 


X-ray bright superbubble DEM 199 


7 


Sk-70°91 


81.89 


-70.61 


14.55 ±0.07 


Superbubble DEM 208 


8 


Sk-66° 100 


81.94 


-66.92 


14.26 ±0.06 


Field star, no Ha 


9 


Sk-67°144 


82.55 


-67.43 


14.41 ±0.07 


Periphery of supergiant shell LMC 4 


10 


Sk-71°45 


82.81 


-71.06 


14.37 ±0.02 


Superbubble DEM 221 


11 


Sk-69°191 


83.58 


-69.75 


14.38 ±0.06 


X-ray bright superbubble DEM 246 


12 


Sk-67°211 


83.80 


-67.55 


14.20 ±0.02 


H II region DEM 241 


13 


BI237 


84.06 


-67.65 


14.39 ±0.03 


Diffuse H II region 


14 


Sk-66° 172 


84.27 


-66.35 


14.31 ±0.06 


H II region DEM 252 


15 


Sk-69°246 


84.72 


-69.03 


14.36 ±0.08 


X-ray bright superbubble 30 Dor 


16 


<N70> 


85.81 


-67.85 


14.48 ±0.10 


Superbubble DEM 301 


17 


CAL83 


85.89 


-68.37 


14.39 ±0.08 


Field star 


18 


BI272 


86.09 


-67.24 


14.30 ±0.03 


Diffuse H II region DEM 309 


19 


Sk-67°266 


86.46 


-67.24 


14.09 ±0.07 


Shell DEM 3 15 



Note: O VI column densities are from this work. lHowk et all fe002l), and Danforth & Blar] 120061) . 
The designations DEM refer to entry in the catalog of lDavies. Elliott. & Meaburnl7l9767i 



components could have been ejected away from the LMC in this 
process. 

- The non - equilib rium radiative cooling model of 
Edga r & Chevalier! |l986) where a gas cools from tempera- 
ture of > 10 6 K fits also the data when O V I/C IV is greater than 
5.9. More recently, iGnat & Sternberg! d2006t) have produced new 
calculations of radiatively cooling gas as a function of metallicity 
using updated atomic data. The O VI/C IV and O Vl/Si IV ratios 
for the HVC and LMC absorbers are consistent with the Gnat 
& Sternberg models of radiative cooling from T ~ 5 x 10 6 K 
for models with metallicities 0.1 to 1.0 times solar. The gas must 
still be at temperatures between ~ (1 — 2) x 10 K in order to 
match their models. A radiative cooling model for the origins of 
the highly-ionised gas in the HVC and LMC components would 
be consistent with an outflow and infall of highly ionised material 
from the disk of the LMC, i.e., an exchange of matter between the 
thin disk of the LMC and its halo or thick disk. The majority of the 
gas must cool in the outflow phase since we only find evidence for 
the possible infall of m aterial towards one sight line. If we combine 
our sample with that o f lHowk et alj J2002I) . only 3 of 14 sight lines 
show possible evidence for infalling material. It is, however, not 
clear within the framework of this model why velocities of the cool 
and hotter gas would be systematically coincident. 

In summary, the physical origins of O VI and broad C IV and 
Si IV are tangled and several processes could occur. However, it 
is apparent that conducting heating processes play a fundamental 
role in the production of the highly ionised species probing gas at 
T ~ 10 5 K while TMLs are not important in the probed regions in 
the HVC-LMC, except if a large fraction of O VI is produced via 
another mechanism (e.g., radiative cooling). 

4.2.4 Fluctuation ofN(OWl) 

We note that lDanforth & Blairl {2006) found that the O VI column 
densities appear higher by about 0. 2 dex in the N70 s uperbubble 
and other superbubbles described bv lHowketaT1 (f2002) compared 



to some more quiescent H II region sight lines, suggesting the very 
local environment of the background stars may affect the total col- 
umn of O VI along a sight line. In this context, they argued that su- 
perbubbles were significant contributors to the total O VI observed 
along LMC sight line s. Fig. 8 shows an Ha image of the LMC from 
iGaustad et al.1 f200lb with a representatio n of all available O VI 
column densities superimposed (including iDanforth & Blairll2006t 
iHowk et al.l2002l , and this work). Table[7]summarises the measured 
columns and the surrounding environments for each FUSE mea- 
surement. All of the superbubble sight lines are concentrated below 
the line superimposed on the image (we display only the average of 
the four measurements towards N70). 

While some superbubble-type regions may have slight en- 
hancements in O VI, what is striking about this figure is the very 
large local variations, as reported by Ho wk et all d2002l) . For ex- 
ample the pairs of stars labeled 12 & 13 (Sk-67°211 & BI 237) 
and 18 & 19 (BI 272 and Sk-67°266) trace H II regions where 
the column densities varies by a factor 2 from 14.1 to 14.4 dex. 
These pairs have projected separations o f 90 and 129 pc, r espec- 
tively. This is the degree of enhancement IDanforth & B lair ( 2006) 
identify with superbubbles, but over much smaller angular scales 
than those authors use for comparison between superbubble and 
non- superbubble sight lines. Furthermore quiescent regions like 
stars no. 2, 5, and 17 have O VI column densities similar or larger 
than observed in some superbubbles. The dispersion of N (O VI) is 
large along superbubble sight lines, as well. The measurements of 
A^(0 VI) appear not to be smaller than 14.3 dex in the superbub- 
bles. However, many more O VI measurements will be needed in 
the LMC before one can truly determine if the superbubbles and 
other very local phenomena have a strong impact on the observed 
total O VI column density. It is clear that the O VI has a patchy 
distribution in the LMC, and comparisons between sight lines sep- 
arated by tens of parsecs may not lead to significant results in this 
regard. While the high-ion ratios vary from sight line to sight line, 
there is no clear link between the observed variation and the regions 
associated with the lines of sight. The comparison of high-ion ratios 
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j • i 1 i i : 

06 h 00 m 05 h 30 m 00 

RA (J2000) 

Figure 8. Ha image of the LMC jGaustad et alj|200ll) . The identification numbers are summarised in Table [7] Darker regions correspond to brighter Ha 
emission. The radius of the circle marking each probe star is linearly proportional to the column density of interstellar O VI at LMC velocities. Note that the 
apparent emission near no. 5 is caused by an imperfectly subtracted foreground star. Fortuitously, the superbubbles are well separated from the H II regions 
(except for no. 17, which is an H II region) and the solid line marks the separation between them. Note that several H II regions are associated with supergiant 
shells. 



is potentially a more secure indicator of the impact of the local ef- 
fects because it is not subject to the varying depth of the gas probed 
by the background stars, which may be part of the source of the 
observed variations in N(0 VI). 



5 SUMMARY 

We analyse the physical properties of the highly ionised species 
(C IV, Si IV, N V, O VI) in and around the LMC towards a sam- 
ple of 4 LMC hot stars. Two of these sight lines probe prominent 
superbubbles, while the other two probe H II regions, allowing us 
to compare the highly ionised species and their connection with 
weakly ionised species among different regions. O VI was obtained 
with FUSE with a spectral resolution of ~20 km s~ 1 , and the other 
highly ionised species were observed with the HST/STIS E140M 
setup with a spectral resolution of ~7 kms -1 . The excellent spec- 
tral resolution and high signal-to-noise of these data have allowed 



us to model the highly ionised species in the LMC. Our main find- 
ings are summarised as follows: 

(i) Absorption from both the high and low ionisation species 
are observed over the entire range of velocities from the Milky Way 
to LMC (-50 < Wlsh < +350 kms -1 ). Along our four sight 
lines, the profiles are separated into four principal components: the 
MW (wlsr < 30 kms" 1 ), IVC (30 < v hS R < 100 kms" 1 ), 
HVC (100 < «lsr < 175 kms" 1 ), and LMC («lsr > 175 
km s~ ). N V is not detected at the 3a level. 

(ii) In the LMC component, we find narrow and broad absorp- 
tion in C IV and Si IV, but only very broad O VI absorption. The 
narrow LMC components of C IV and Si IV have an observed b- 
value that implies temperatures of a few times 10 4 K or less. The 
properties of the narrow components of C IV and Si IV can be prin- 
cipally explained by the interstellar environment associated with 
the stars, in particular H II regions around the central star and ex- 
panding shells behind the shocked plasma. 

(iii) The O VI at LMC velocities is broader than the O VI 
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absorption observed in the combined Milky Way and IVC compo- 
nent, implying that the LMC O VI profiles are composed of several 
components. The broad LMC components of O VI, C IV, and Si IV 
likely probe gas of a few times 10 K. 

(iv) We find a striking similarity in the O VI/C IV ratios for the 
LMC and HVC components (wlsr > 100 km s~ ) along a given 
sight line, but the O VI/C IV ratios for the HVC and LMC gas 
vary between the sight lines (see Fig. 7). The fact that the ratios in 
the HVC and LMC components along a given sight line are always 
very similar suggests that the HVC component has its origins in 
the LMC, and the processes producing the highly-ionised gas in 
the LMC affect a large range of velocities along a given sight line. 

(v) The difference in the high-ion ratios between the four sight 
lines implies different processes or varying conditions may produce 
the highly ionised species in the LMC. There is no evidence, how- 
ever, of a link between the observed variation and the regions as- 
sociated with the stars, although it is likely that some of the O VI 
absorption occurs because of stellar winds or other phenomena as- 
sociated with the stellar environment. Conductive interface model 
can reproduce both the high-ion ratios observed in the broad com- 
ponent and the apparent kinematically coupling between O VI and 
the weakly ionised species such as Si II and O I. 

(vi) The velocities of the LMC-IVC and HVC components in- 
dicate that they are expanding away from the warm disk gas. Only 
towards one line of sight is there evidence for both significant red- 
shifted and blue-shifted components in the low- and high-ion spec- 
tra, implying little of the material we observe is returning to the 
LMC disk. Our analysis therefore provides compelling evidence 
for a hot LMC halo fed by energetic outflows from the LMC disk 
and even possibly with a galactic wind since the velocity of the 
HVC relative to the LMC disk is actually large enough to escape 
the LMC. 
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